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Abstract 

The possibility of realizing the superradiant regime of electromagnetic emission 
by the assembly of quantum dots is considered. The overall dynamical process is an- 
alyzed in detail. It is shown that there can occur several qualitatively different stages 
of evolution. The process starts with dipolar waves triggering the spontaneous radia- 
tion of individual dots. This corresponds to the fluctuation stage, when the dots are 
not yet noticeably correlated with each other. The second is the quantum stage, when 
the dot interactions through the common radiation field become more important, but 
the coherence is not yet developed. The third is the coherent stage, when the dots 
radiate coherently, emitting a superradiant pulse. After the superradiant pulse, the 
system of dots relaxes to an incoherent state in the relaxation stage. If there is no 
external permanent pumping, or the effective dot interactions are weak, the system 
tends to a stationary state during the last stationary stage, when coherence dies out 
to a low, practically negligible, level. In the case of permanent pumping, there ex- 
ists the sixth stage of pulsing superradiance, when the system of dots emits separate 
coherent pulses. 



PACS numbers: 73.21.La, 73.21. Fg, 78.67.Hc, 78.67.De 



1 Introduction 



Superradiance is the effect of self-organized collective coherent radiation by an ensemble of 
radiators. The phenomenon of optical superradiance is well known for atoms and molecules, 
being described in numerous publications (see., e.g., books [1-3]). Optical superradiance 
from Bose-Einstein condensed atoms has also been investigated [4]. There also exists 
spin superradiance produced by spin systems, such as nuclei [5-8] or magnetic molecules 
[9-12]. It has been suggested [13] that the assemblies of quantum dots or wells could 
also be arranged so that to produce superradiance. This can become possible when the 
distance between the neighboring quantum nanostructures is smaller than the radiation 
wavelength. In that case, there appears an effective interaction between the radiators, due 
to common radiation field. The modern technology of preparing materials with quantum 
dots allows for the fabrication of the quantum dot assemblies with the density of quantum 
dots sufficient for the appearance of such an effective interaction [13-17], which has been 
detected experimentally [18]. 

Quantum dots, whose electrons or holes are confined in all three spatial dimensions, 
have many properties [19-22] that make them similar to atoms, because of which quantum 
dots are often called "artificial atoms". One of the main such properties is the existence 
of discrete energy levels, whose shell structure can be adjusted by design. Some of these 
properties are shared by quantum wells confining electrons or holes in one dimension and 
allowing free propagation in two dimensions. For concreteness, we shall concentrate in 
what follows on the consideration of quantum dots. 

Transferring electrons from the ground-state energy level to an excited level creates a 
hole. The interacting pair of an excited electron and a hole forms an exciton. The electron- 
hole recombination is accompanied by the radiation of electromagnetic field, which is in a 
very close analogy with the radiation of excited atoms. This is why it was reasonable to 
assume [23] that an ensemble of quantum dots can be employed for the creation of quantum- 
dot lasers. Several types of quantum-dot lasers have been demonstrated since then, based 
on semiconductor heterojunctions [24-29] and photonic crystals [30-33]. Quantum- well 
lasers have also been realized. But the latter, because of thermal occupation of the the 
quasi-2D continuum, are essentially more sensitive to temperature changes, which makes 
their emission less monochromatic than that of quantum dots. Quantum-dot lasers also 
enjoy a lower current- density lasing threshold. Because of these differences, quantum dots 
look more suitable for the use as radiating devices. 

When quantum dots are fabricated in the process of epitaxial growth, they are usually 
characterized by a noticeable size dispersion, which results in the related inhomogeneous 
broadening. The latter may hinder the possibility of achieving efficient dot interactions 
through the common radiation field and, as a result, destroying the cooperative character of 
emission, thus, suppressing coherence [15,34,35]. However the impressive recent progress in 
controlling quantum-dot parameters [18,27] makes it now feasible to create quantum dots 
so that they allow for the appearance of effective dot interactions through electromagnetic 
field, which is a prerequisite for achieving collective coherent radiation. 

Because of the technological feasibility of fabricating semiconductor samples with suffi- 
ciently dense and uniform quantum dots, it should be possible to realize the conditions for 
their superradiant emission. It is the aim of the present paper to consider the quantum- 
dot superradiance. The main goal is to develop an accurate and detailed description of all 
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stages of the superradiant dynamics, starting from the initial stage, when collective effects 
are yet weak, to the coherent stage of radiation, and further to the end of the whole pro- 
cess. The developed theory makes it possible to describe different types of superradiance, 
such as pure superradiance, triggered superradiance, and pulsing superradiance. 

The aim of the paper is to consider the realistic situation corresponding to quantum 
dots, but not an oversimplified model consideration. Therefore, in order to understand 
what approximations are admissible for characterizing the process, in the next section, an 
analysis of the parameters, typical of semiconductors with quantum dots, will be given. 
Such an analysis is necessary before plunging into mathematical formulas in the following 
sections. The parameters are taken from the literature cited above. 

2 Typical Material Parameters 

For realizing quantum-dot superradiance, it is reasonable to take those materials that 
are used for quantum-dot lasers. For the latter, one usually employs the self-assembled 
heterostructures, such as InAs/GaAs, InGaAs/GaAs, InGaAs/AlGaAs, GalnAsP/InP, 
InAs/InP, InAs/GalnAs, AlInAs/AlGaAs, InP/GalnP, AlGaAs/GaAs, and CdSe/ZnSe. 
There exist different kinds of quantum dots, having different shapes and sizes. Thus, the 
lateral size of a typical self-assembled quantum dot is much larger than its vertical ex- 
tent. Typical dot sizes are of the order r dot ~ 10~ 7 cm — 10~ 6 cm. In each dot, there can 
be between just a few to 10 5 electrons. The dot density in an epitaxy layer is of order 
10 8 cm~ 2 — 10 n cm -2 . With the width of the layer h ~ 10 _6 cm — 10 _5 cm, this makes the 
spatial density p ~ 10 13 cm~ 3 — 10 17 cm -3 . The interdot distance is a ~ lCT 5 cm — 10~ 4 cm. 
The lasing operation is realized at the wavelength A ~ 10~ 4 cm, which translates into the 
frequency ooq ~ 10 15 Hz. The natural width 70 = 2|g?| 2 A;q/3, where ko = 2ir/X, depends on 
the transition dipole d. The typical value for the latter is d ~ 100D. Taking into account 
that ID = lO^^/en? • cm 3 gives 70 ~ 10 10 Hz. The actual homogeneous broadening is 
usually larger than 70, being 72 ~ 10 12 Hz — 10 n Hz. For high-quality self-assembled dot 
materials, the inhomogeneous broadening can be made relatively small, of the same order 
as 72, that is, of order 72 ~ 10 12 Hz — 10 13 Hz. The longitudinal relaxation time 7i is mainly 
due to electron-phonon coupling, which can be strongly suppressed by low temperatures. 
Thus, at helium temperatures, the longitudinal dephasing time becomes limited only by 
the lifetime of inversion for a single quantum dot in free space, Ti ~ lCT 9 s, which gives 
7! = 1/Ti ~ 10 9 Hz. To enhance a chosen mode, one places the sample into a resonator 
cavity with a large quality factor reaching 10 4 . The sizes of the sample can be different. 
Edge-emitting lasers do not have a circular cross-section. The typical sizes of quantum dot 
lasers can be R ~ 10~ 3 cm — 10~ 2 cm. 

Since the radiation wavelength is much larger than the dot sizes, A ^> r dot , the dipole 
approximation is appropriate. For sufficiently dense dot materials, the interdot distance 
can be made much smaller than the wavelength, a < A. Hence, there can arise sufficiently 
strong dot interaction through the common radiation field. But the wavelength is much 
smaller than the sample linear sizes, A < I. Therefore the Dicke-type [36] approximation 
of a concentrated, point-like, sample, cannot be used. In strongly nonuniform materials, 
with a very large inhomogeneous broadening 7I, such that T 2 * = I/72 is comparable with 
the time of the radiation pulse t p , superradiance is suppressed [34]. It is therefore neces- 
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sary to prepare the samples with a narrow distribution of dot sizes, resulting in not too 
wide inhomogeneous broadening. Fortunately, the fabrication of such samples is nowadays 
technologically possible. As is seen from the above values for typical parameters, the in- 
homogeneous broadening can be made of the same order as the homogeneous one. When 
72 ~ 72, then the consideration can be simplified by combining ^ an d 72 m one effective pa- 
rameter [37]. Thus, the typical dephasing time is rather short, T 2 = I/72 ~ lCT 13 s— lCT 12 s. 
Superradiance is possible only if the time of radiation pulse r p is shorter than the dephasing 
time T2. In the following sections, the theoretical description is developed, which takes into 
account the typical characteristics of the quantum dot assemblies. Throughout the paper, 
the system of units is used, where the Planck constant h = 1. 



3 Basic Operator Equations 



Aiming at developing a realistic description of the system, let us start with the microscopic 
Hamiltonian 

H = H d + Hj + H d j + H m f , (1) 

characterizing an ensemble of radiating quantum dots in a semiconductor matrix inside a 
resonator cavity. The Hamiltonian 

N 

(2) 



where Sf is a pseudospin operator of an i-th dot, represents N two-level quantum dots, 
with the carrying transition frequency uq. Mathematically, the operator Sf is just a spin 
operator. It is called the pseudospin operator, since it corresponds not to an actual spin 
but to the population difference. The radiation-field Hamiltonian 

H f = ^J (E 2 + H 2 ) dr, (3) 

contains electric field E and magnetic field H. The vector potential A, introduced by the 
standard relation H = V x A, is assumed to satisfy the Coulomb calibration V • A = 0. 
The dot-field interaction is given by the Hamiltonian 



N (\ \ 
II -]T -J,- A, • IVKo, • (4) 



l df 



in which Aj = A(r i: t) and E 0i = ~E (r i: t) is a seed electric field of the resonator cavity. 
The transition current 

J, = iu (dS+ - d*Sr) , (5) 

and transition polarization 

Pi = dS+ + d*S~ , (6) 

are expressed through the transition dipole d and the ladder operators Sf = Sf±iSf : where 
Sf = S a (ri,t). Since the cavity is filled by a semiconducting material, the Hamiltonian 

H m f = - ~ c J jmat(r, t) ■ A(r, t) dr , (7) 
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describes the interaction of the local density current j ma t( r , t) in the filling matter with the 
radiated electromagnetic field. 

The field operators satisfy the equal-time commutation relations 



[E a (r,t), AP(r',t)} = 4mc5^(r - r') , 

d 



[E a (r,t), HP(r',t)} = -Amc^a/h ^~ 6(r - r') , (8) 
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where e Q( g 7 is the unitary antisymmetric tensor [3] and the transverse delta-function is 

Mr) - / (<W " ^f) e lkr ^ = I MOO - -L D a p (r) , (9) 
with the dipolar tensor 

_ , x 8 a p - ?>n OL n fi 
D a p{r) = - 



in which n = r/r = {n a } and r = |r|. The pseudospin operators satisfy the spin commu- 
tation relations 

[5* Sr] = 2^-5? , [5?, Sf\ = ±5 l3 Sf . (10) 
The Heisenberg equations of motion for the field operators yield 
1 <9E 47r . 1 <9A 



with the density of current 



N 

tapir - Ti)jf(t) + / 6 a p(r - r')/ mat (v', t) dv> 
i=i ^ 



(12) 



^From these, using the Coulomb calibration, one gets the equation for the vector potential 

The Heisenberg equations for the pseudospin operators result in the equations 
dS- 



dt 

dSf 



= -tuj Sr + 2S* (k d • Aj - id • E w ) , 



-5+ (A; d • Aj — id • E 0i ) - S~ (k d* • A, + *d* • E w ) , (14) 

where ko = ojq/c. These equations are to be complimented by the retardation condition 

S?(t) =0 (t<0). (15) 

Equations (11) to (15) are the basic operator equations describing all radiation processes 
in the system of quantum dots. It is worth emphasizing that these equations follow from 
the first principles. Such a microscopic approach is necessary for correctly treating the 
overall dynamics of quantum dot radiation. 



5 



4 Elimination of Field Variables 



The standard way of considering the radiation processes is by averaging Eqs. (11) and (14) 
and passing to the semiclassical approximation. This way, presupposing well organized 
coherence, does not allow for the treatment of those radiation stages, when coherence has 
not been developed. Therefore, we employ here another, more accurate, approach allowing 
for the treatment of all radiation stages. 

The first step in the approach to be pursued is the elimination of field variables [8]. 
This can be done by solving Eq. (13) for the vector potential and substituting the found 
solution into Eqs. (14) for the pseudospin operators. The known solution of Eq. (13) is 
the sum 

A(r, t) = A vac + - J j [r>, t - j-—^ (16) 

of the vacuum potential and the retarded potential, with the density of current given by 
Eq. (12). This solution is to be substituted into Eqs. (14). 

The interaction between the radiation field and a radiating dot is assumed to be small 
as compared to the transition frequency. This is a necessary requirement for the existence 
of well defined energy levels in a dot. In the other case, the transition frequency would not 
be defined in principle. This implies that the retardation in the time dependence of the 
current density can be taken into account in the Born approximation 

ST (t-^= S-(t)Q(ct - r)e ikor , S* (t - = S*(t)Q(ct - r) , 

where 6(t) is the unit step function. 

When substituting the vector potential (16) into Eqs. (14), one meets the terms cor- 
responding to the dot self-action. The contribution of these terms is characterized in the 
Appendix. Finally, the vector potential (16) is represented as the sum 

A = A mc + A. se if + A mc ; + A^p + A ma £ . (17) 

Here the first term is due to vacuum fluctuations. The self-action potential is described in 
the Appendix. The radiation potential 

A ^.') = E^^ J ,( i -^) • < 18 > 

in which r ^ rj, describes the spherical part of the vector potential, corresponding to the 
radiation field produced by dots. In addition, the radiating dots create the dipolar part of 
the vector potential 

W^-££/%^4-^)*<. (19, 

The interaction of the semiconductor, filling the cavity, with the radiation field produces 
the potential 

4Ur, *)=£/ ^TZ^P ?™* ( r "> * " ^) dr ' dr " ■ ( 20 ) 
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Substituting the vector potential (17) into Eqs. (14), we obtain the equations for the 
pseudospin operators, which do not contain the field variables. Instead, there appear 
effective dot interactions through the common radiation field. 

5 Stochastic Mean-Field Approximation 

The dynamics of the system can be characterized by the behavior of the following functions. 

The transition function 

u(r i ,t) = 2{S-(t)) H (21) 

where the angle brackets imply quantum statistical averaging associated with the system 
Hamiltonian H, describes the local effective polarization corresponding to dipole transi- 
tions. The coherence intensity 

9 N 

W ^ = N^2 [(SmSi(t)) H + (Sf(t)Sr(t)) H ] (22) 

is the local characteristic of coherence. The intensity of coherent radiation is proportional 
to this function. The population difference 

s(r h t) = 2(Sm)H (23) 

defines the local difference of populations for the energy levels of an i-dot. 

To obtain the evolution equations for functions (21), (22), and (23), we have to average 
the equations resulting after the substitution of the vector potential (17) into the equations 
of motion (14). Such equations are not closed. To make them closed, it is necessary 
to invoke some decoupling for the operator correlation functions. If one resorts to the 
standard mean-field decoupling, one comes to the usual semiclassical approximation. As 
is well known, the semiclassical approximation can be used only when the system is in a 
coherent state or is almost coherent [2,3]. But incoherent regimes cannot be described in 
this approximation. One of the most interesting questions is how coherence develops in an 
initially incoherent system. To be able to describe such a regime, it is necessary to employ 
a more accurate approximation. For this purpose, we shall use the stochastic mean-field 
approximation employed earlier for describing the dynamics of spin assemblies [5-11] and 
Bose systems in random fields [38-40]. 

Let us combine the vacuum, dipole, and matter vector potentials into the sum 

£(r, f) = 2k d ■ (A mc + A dip + A mat ) . (24) 

The potentials, entering this sum, create local fluctuations of electromagnetic field. Being 
averaged over space, quantity (24) is practically zero. Therefore Eq. (24), characterizing 
the strength of local field fluctuations, can be treated as a local random variable. Contrary 
to this random variable, the radiation potential (18) induces long-range effective interac- 
tions between the radiating dots. The long-range radiation potential (18), with current 
(5), is expressed through the pseudospin operators Sj. Hence the effective interactions, 
arising between dots, are mathematically similar to long-range spin interactions. 
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There is a direct similarity between the effective dipole interactions caused by the dipole 
vector potential (19) and the dipole [5-11] and hyperfine interactions [9-11,41,42] in spin 
systems, these interactions being treated as stochastic fluctuations playing destructive role 
by dephasing collective motion. While the effective pseudospin interactions, induced by 
the radiation potential (18) in quantum dots, are equivalent to the effective interactions 
produced by the resonator feedback field in spin systems [9-11]. The latter interactions are 
responsible for the appearance of collective effects in magnetic and ferroelectric samples, 
including the arising coherence [43]. In spin systems without a resonator feedback field, 
coherence and, hence, superradiance, cannot develop, being destroyed by the direct dipole 
interactions [5-11,44,45]. 

In this way, it is possible to distinguish two types of variables in the system, the random 
variable £ describing local field fluctuations and the pseudospin variables Sj characterizing 
effective long-range dot interactions. For brevity, we may denote the collection of pseu- 
dospins {Sj : j = 1,2, ...,7V} by S. Then, the operators of observable quantities O are, 
generally, functions of these variables, O = 0(S,£). 

Having two kinds of variables, it is possible to define two types of averaging. One type 
corresponds to the quantum statistical average 

(d) H = Trp6(S,0 (25) 

involving the pseudospin operators S, with p being the system statistical operator, and 
where the random variable £ is kept fixed. Another type is the stochastic averaging over 
the random fluctuations, denoted as 

«0» = Jp6(8,t)Vt, (26) 

which is defined through a functional integral over the random variable £, with the pre- 
scribed differential measure V£. Respectively, the total averaging 

(0) = {{{{0) H ))) (27) 

includes both, the quantum and stochastic, averages. 

Keeping in mind the long range of the effective pseudospin interactions, we can use the 
stochastic mean- field decoupling 

(S?Sf) H = (S?) H (Sf) H (i^j), (28) 

where only the quantum averaging is involved. This decoupling looks like a mean-field 
approximation. However, it has a very important principal difference form the latter, 
involving only the quantum averaging (25), but not touching the stochastic averaging (26). 
Since no approximation is done here with respect to the stochastic variables, decoupling 
(28) preserves stochastic properties of the system. This is why it is called the stochastic 
mean-field approximation [5-11]. 

When the radiation wavelength A is much larger than the interdot distance, the geo- 
metrical location of dots in space is of no importance. Then it is convenient to pass to the 
continuous spatial representation, replacing the sums by the integrals according to the rule 

j^^pjdv (p=y), (29) 
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with the integration over the whole system and p being the dot density. 

In order to represent the evolution equations in a compact form, let us introduce the 
effective field, or effective force, acting on dots 

/(r, t) = / (r, t) + f rad (r, t) + £(r, t) . (30) 

Here, the first term 

/o(r,f) = -2id.E (r,*) (31) 
is due to the external field E . The second term 

f rad (r,t) = 2k Q (d-A rad (r,t)) H (32) 

is caused by the radiating dots. And the last term in Eq. (30) is the fluctuating random 
field (24). The radiation term (32), with the vector potential (18), acquires the form 



frad(r,t) = -i-fop J 



d 



2 



G(r - r', t)u(r', t) - G*(r - r', t)u*(r', t) 



dr' , (33) 



in which the transfer kernel is 



G( M )E^e(d-r). (34) 

Finally, we average Eqs. (14) according to the quantum averaging (25), employ the 
stochastic mean-field decoupling (28), and use notation (30). Then for the variables (21), 
(22), and (23), we obtain the evolution equations 

— = -{iu + 7 2 )m + js, — = -2^ 2 w + {u f + / u)s , 
(9s 1 

^ = -2(«V + r«)-7i(«-0, (35) 

in which u = u(r,t),w = w(r,t), and s = s(r,t). The longitudinal, 71, and transverse, 72, 
attenuation rates are treated as the system parameters, whose typical values are discussed 
in Sec. II. The parameter ( characterizes the level of stationary nonresonant pumping. 
The evolution equation for w follows from definition (22), with the use of the Heisenberg 
equations of motion for the pseudospin operators and the stochastic mean-field decoupling 
(28). Since in definition (22), the summation is over j 7^ i, it is also possible, first, to 
decouple the products of the pseudospin operators, according to Eq. (28), and then invoke 
the Heisenberg equations. In any case, the result is the same equation for w. 

Equations (35) are stochastic integro-differential equations. These are the basic evolu- 
tion equations describing the dynamics of radiation in a system of quantum dots. 



6 Triggering Dipolar Waves 

As is mentioned above, the dipolar vector potential (19) induces local fluctuations dephas- 
ing the radiation of the dot assembly. However, these fluctuations play not only destructive 
role. They can be useful at the initial stage of the radiation process, when the latter is not 
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triggered by an external field. In such a case, the transition dipole fluctuations, related to 
spontaneous emission, can trigger the process of collective radiation. 

In order to illustrate the appearance and the nature of the dipolar fluctuations, let us 
consider the pseudospin equations (14), leaving there only the terms related to the dipole 
interactions. Let us introduce the interaction coefficients 

.2 J.2 



a/3 a/3 



where 



D°f = / Q(ct - \ Ti - r'|) Daf){r '_ rj) exp(-^ |r t - r'|) dr> . 



The latter quantity, because of the unit-step function 0(ct — |r» — r'|) in the integral, varies 
in time only at the very beginning of the process, when i < o/c, after which it becomes 
practically constant. Therefore, for all times, except f C a/c, the interaction coefficients 
(36) can be treated as constant parameters. 

For the dipolar vector potential (19), we have 

m • A dip ( ri ,t) = -^Yl M/w - ^nt)] ■ 

Then the equations of motion reduce to 
dS7 



dt 

Hi*) 

~£ = 2 £ l S t M - C - S i) - S i {^Sj - 45;)] . (37) 
The dipolar pseudospin fluctuations are described by the small deviations 

8Sf = St- (Sf) , 8S? ee St - {St) (38) 

from the related average values < Sf >. The latter averages correspond to equilibrium or 
quasiequilibrium, when they either do not depend on time or are slow functions of time, 
as compared to the fastly fluctuating deviations (38). Linearizing Eqs. (37) with respect 
to small deviations (38), under the condition < S~ >= 0, and taking into account that, 
because of the properties of dipole interactions, 

£ K = £ % = , (39) 

we obtain the equations 

d 5Sr = -iu 5Sr - i{St) J2 (MS/ - CijSST) , 



dt 



|^ = 0. (40) 
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Then we employ the Fourier transforms for the pseudospin operators 

S t = E St ex P(^k • r,-) , Sf = 1 £ Sf exp(±*k • r,) , (41) 



AT 



for the interaction coefficients 

1 X 

% = E &fc eX P( ik ■ > ^ = E 6 « eX P( _ik ' » ( 42 ) 

and, similarly, for the interaction coefficients Cy, where = — r^. Taking into account 
that 

SSf = Sf ((Sf) = 0) (43) 

and introducing the notation 

^ = w - Cfc(S'f) , A fe = 6_ fc (^), (44) 

we come to the equations 

^ = -vi fc S fc " - iA fc S+ , -JL = ifiS+ + iKS; . (45) 

Looking for the solutions to Eqs. (45) in the form 

S k = u k e~ lu)kt + v* k e iu>kt , S+ = u* k e iuikt + v k e~^ kt , (46) 
we find the dipolar-wave spectrum 

Wfc = \/KI 2 - \M 2 ■ (47) 

This means that the dipolar part (19) of the vector potential generates local field fluctu- 
ations, having the meaning of the transition dipolar waves with the spectrum (47). The 
interaction coefficients (36) are smaller than the transition frequency, so that 

(48) 

Hence, 

IX, I IX, I 

(49) 

Therefore the spectrum (47) is always positive, implying that the dipolar waves are dy- 
namically stable [46]. In the long-wave limit, when k — >• and 

1 y 1 . 

hk - ~~ 2 ^ bij(yk ' ^ ' Ck ~ ~ 2 ^ Cij(yk ' ' 

j(#0 J'(#0 

the spectrum becomes quadratic, which follows from the expression 

■ 2 ~ w 2 + u (S?) Re J2 c u( k ' r u) 2 • (50) 





«1, 




< 1 






w 


w 






«1, 




< 1 






w 







The existence of these dipolar waves triggers the process of dot radiation, even if no external 
field is imposed at the initial time. 



11 



7 Transverse Mode Expansion 



After the radiation process is triggered by the dipolar waves, the overall radiation dynamics 
is described by Eqs. (35) that are stochastic integro-differential equations in partial deriva- 
tives. If the sizes of the whole sample would be much smaller than the radiation wavelength, 
we could essentially simplify the problem by resorting to the concentrated-sample approx- 
imation [1-3,36], when just one mode fills the cavity. But in realistic situation, vice versa, 
the wavelength is usually much smaller than the sample sizes, so that the sample can house 
several modes. Then, to simplify the equations, it is necessary to specify the shape of the 
resonator cavity. 

For concreteness, let us assume that the cavity is cylindrical, with radius R and length 
L, such that the wavelength is much smaller than these sizes: 

A « 1 , A « 1 . (51) 

Directing the cylinder axis along the axis z, we can treat the resonator seed field as prop- 
agating along this axis, that is, having the form 

E (t, t) = - E x e i(kz - Wt) + - £* e -^-"') . (52) 

The cavity is resonant in the sense of the small detuning of the resonator natural frequency 
from the dot transition frequency, 

^ < 1 (A = oj — ojq) . (53) 

COq 

In the same way as it is for all dot-field interactions, for the seed field, we have 



l^il 
u 



<1 (i/i = |d-Ei|) . (54) 



In the cylindrical geometry, the radiation modes acquire the shape of filaments extended 
along the axis z. Then, it is possible to represent the solutions to Eqs. (35) as expansions 
over the transverse modes: 

N f N f 

u(r, t) = Y^ u ni r ±, t)e ikz , w(r, t) = ^ w n (r±, t) , 

n=l n=l 



N 



/ 



S M) = ^2 s n(r±,t) , (55) 



n=l 



where Nf is the number of the filamentary modes and r± = a/x 2 + y 2 is the transverse 
radial variable. Ascribing to a mode an effective enveloping radius Rf, for the effective 
enveloping volume of a filamentary mode, we have Vf = nRjL. Representation (55) is 
rather general, including the case of just a single mode, when Nf = 1 and Vf = V, with V 
being the sample volume. We keep in mind that the sample and cavity volumes coincide. 
This assumption does not reduce the generality of consideration, since when these volumes 
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are different it is sufficient to take into account the existence of a filling factor that is not 
equal to one. 

In what follows, we assume that the distance df between the axes of any two nearest- 
neighbor filaments is larger than twice the filament enveloping radius. The physical condi- 
tion, corresponding to this assumption can be understood by taking into account that the 
filament enveloping radius is of the order of the diffraction radius, i.e., Rf ~ y/\L. Then, 
from the inequalities Rf <^ df < R, it follows that the Fresnel number F = R 2 /\L is to 
be large, F 1. Therefore, in the case of large Fresnel numbers, different filamentary 
modes can be treated as uncorrelated and considered separately from each other. In such 
a case, we can reduce the consideration to studying the behavior of each mode on average 
by defining the averages over the enveloping volume of a mode as 



If 2 f f 

u(t) = — / Un(r±,t) dr = u n (r,t)rdr , 

Vf Jvt Rf Jo 

If 2 f Rf 

M7 ( t ) = 77" I w n( r ±,t) dr =~B2 w n (r,t)rdr , 
v f JVf K f Jo 

If 2 f R f 

s(t) = — / s n (r±,t) dr= —2 s n (r,t)rdr , (56) 
v f Jv f K f Jo 

where, for the simplicity of the following notation, the index n enumerating modes is 
omitted in the left-hand side of these equations. 
We also need to introduce the coupling functions 



Top / e(ci-r) E^zM rfr , 

Jv f k r 
P{t) - 70P J Q{ct - r) C ° S( ^ ~ kz) dr , (57) 



'v f 

the average stochastic field 

m = yj e(r,t)e^ 2 rfr, (58) 

and the effective force 

= -kl -E ie --* + £(*)■ (59) 

Then, we substitute the mode expansions (55) into Eqs. (35), averaging the mode functions 
according to Eqs. (56). We use the condition that different spatial modes are not correlated 
with each other, so that 

J2 u m(r±,t)e ikz s n (r ± ,t) = ^ u n (r±, t)e ikz s n (r ± , t) . 

mn n 

And we employ the theorem of average in the integral 

/ G(r' - r, t)u n (r ± , t)e tkz s n (r' ± , t) drdr' = u(t)s(t) [ G{r' - r, t)e ikz drdr' . 

JVf JVf 
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Thus, using the above notations (57), (58), and (59), we obtain the equations 

du 

— = -i{co + (3s)u - (72 - as)u + f ± s , 

^ = -2(72 - as)w + (u*/i + flu) s , 
ds 1 

- = -aw- - (u*h + flu) s - 7 i(s - C) , (60) 

describing the evolution of an averaged mode characterized by functions (56). 

In this way, the mode expansions (55) allow us to transform Eqs. (35) in partial 
derivatives to Eq. (60) in ordinary derivatives. The used mode-expansion method is based 
on the idea of the eikonal approximation [47,48] . 

Formally, the used mode-expansion method reduces the general problem to a collection 
of effective single-mode radiation problems, each decoupled from the others. Such a re- 
duction is, actually, the main idea of the eikonal approximation. Greatly simplifying the 
consideration, this reduction leaves aside the question of what would be the distribution 
of filament sizes. The latter is a separate problem, depending on the sample shape. 



8 Scale Separation Approach 

Equations (60) are stochastic differential equations that are not easy to solve. Fortunately, 
they can be further simplified using the scale separation approach [5-8,43] that is a variant 
of the averaging technique [49]. In the present section, we employ this approach [5-8,43]. 

It is possible to notice that there are different time scales in this system of equations. 
The attenuation rates, as discussed in Sec. II, are small as compared to the transition 
frequency, thus, defining the small parameters 

^<1, ^-<1, ^<1. (61) 

UJq COq COq 

It, therefore, follows from Eqs. (60) that the function u(t) is fast in time, as compared to 
the slow functions in time w(t) and s(t). It is convenient to introduce other slow functions 
having the meaning of the collective width 

r = 72 - as , (62) 

collective frequency 

D, = uo + /3s , (63) 

and effective detuning 

5 = tu-n = A-l3s. (64) 

These slow functions play the role of quasi-integrals of motion, or quasi-invariants, for the 
fast function u(t). The first of Eqs. (60) can be solved by keeping fixed the quasi-integrals 
of motion, which gives 



A(Oe- (<n+r)(t-t,) ■ (65) 
Jo 
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The seed field (52) is defined up to a phase factor. Then, without the loss of generality, 
the global phase of the seed field can be chosen such that the value u d ■ ~Ei be real, where 
u = u(0). 

The found solution (65) has to be substituted into the second and third of Eqs. (60) 
for the slow functions w(t) and s(t). The right-hand sides of the latter equations are to 
be averaged over the explicitly entering time and over the stochastic variable according to 
the rule 

hm - f «. . .))dt , (66) 

with keeping fixed the quasi-invariants. The stochastic variable £, describing local field 
fluctuations, by its definition, is zero-centered, such that 

«£(*)» = . (67) 

And the correlation function <C £*(t)£(tf) ^> defines the dynamic attenuation rate 

73 = Relim - T dt f ((C{tW)))e^ n+r ^-^ dt' , (68) 

caused by these random fluctuations. The dynamic attenuation rate (68) is essentially 
defined by the currents in the semiconductor sample. These currents are usually much 
stronger than those fluctuating in free space. This fact makes the principal difference 
between the considered case of quantum dots in semiconductor and atoms in free space. For 
the latter, the attenuation rate (68) is usually much smaller than the transverse relaxation 
rate 72, while for semiconductors, on the contrary, 73 > 72- 
Let us introduce the effective attenuation rate 

Ts = 73 + p^r 2 (1 - e~ n ) , (69) 

where \5\ < \T\. Following the described averaging procedure, we obtain the equations for 
the guiding centers 

dw ds 

— = -2(72 - as)w + 2T 3 s 2 , _ = - au ,-r 3 s-7i(s-C). (70) 

Thus, the fast variables have been averaged out, while the derived Eqs. (70) characterize 
the evolution of the slow variables. 



9 Dynamics of Dot Radiation 

The temporal evolution of dot radiation through transverse modes is described by Eqs. 
(70). The solutions to these equations essentially depend on the behavior of the coupling 
functions (57), whose values vary with time. It is possible to distinguish several qualita- 
tively different stages of evolution. 

A. Fluctuation Stage 



15 



At the initial time t — 0, the coupling functions (57) are zero. They remain small 
during the time interval 

< t < t int , (71) 

before the interaction time ti nt = a/c, when the dots have not yet been correlated by means 
of the photon exchange. The interaction time, for the interdot distance a ~ 10~ 5 cm — 
lCT 4 cm is ti n t ~ 10 _15 s — lCT 14 s. In the time interval (71), the radiation process starts, 
being triggered by dipolar waves considered in Sec. VI. These waves correspond to the 
random local field fluctuations. As is seen from the above estimates, the fluctuation stage 
is rather short. During this stage, the functions w(t) and s(t) do not essentially change, 
so that w(t int ) ~ w(0) and s(t int ) s(0). 

B. Quantum Stage 

The quantum stage comes after the time ti nt , when the dot interactions through photon 
exchange come into play, but dots are not yet sufficiently correlated in order that coherence 
would develop. At this incoherent stage, dots radiate independently. The stage lasts till 
the coherence time t co h that is necessary for developing coherence. So, the temporal interval 
related to the quantum stage is 

tint < t < t coh . (72) 

During the interaction time t int , the coupling functions (57) quickly grow, reaching, after 
ti n t, their maximal values: 

<x(t) -> £72 , P(t) ->■ £72 (* > Unt) ■ 

Here we have introduced the dimensionless coupling parameters 

7o f sm(k r - kz) 

9 = P— T dr 7 3 

72 Jv f k or 



and, respectively, 



7o f cos(/c r - kz) 
g = p— / dr. 74 

72 Jv f 



At this stage, the collective width (62) becomes 

r = 72 (l-<7s) (75) 

and the collective frequency (63) is 

Q = u + £7 2 s . (76) 

The role of the resonator seed field (52) is to select the resonant frequency, but its amplitude 
is small, such that |ztJ ^72- Therefore the effective attenuation (69) simplifies to T 3 ~ 73. 

Using the above expressions and considering the case when at the initial time no co- 
herence is imposed by external fields, so that w = w(0) = 0, from Eqs. (70) we have the 
equations 

~dt = 2l3s2 ' dt = + 73 ^ S + ' ^ 
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The second of these equations yields the population difference 



(so ~ -^-) exp{-( 7l + 73 )t} + -J^- . (78) 
V 71 + 73/ 7i + 73 



At short times, when (71 + 73 )£ <C 1, the population difference is 

7lC A 1 _, \j. 1 1 ( _ 7iC \ . ,2,2 



s ~ s - so (71 + 7s)f + o so (71 + 7s) r . (79) 

V 7i + 73/ 2 V 71 +73/ 

Then the coherence intensity behaves as 

o-. „ 1 \ Z' „ 7lC \ 



iu ~ 27 3 s t - 27 3 s (7i + 73) so : t + 

71 + 73/ 

2 „, („, , „, \2 ( „ 7iC "\ f n „ 7iC "\ .3 



+ 77 73(71 + 73) I So ~~ j 2s --^- )?. (80) 
3 V 7i + 73/ V 71 + 73/ 

/,From Eq. (79), it follows that the population difference does not essentially vary during 
this stage, being close to so- Taking into account that 71 <C 73, we find that the coherence 
function (80) is either linear in time or cubic in time, 

w ~ 2 73 sgt (s ^ 0) , 

w ^ ^ihsC't 3 (s = 0) , (81) 

depending on whether there exists or not the initial polarization s = s(0). The value 
of the function w in the former case, if s ~ 1, is much larger than the value of w in 
the second case. In the later case, the value of w is practically negligible. This shows 
that in order that coherence could really develop, it is necessary to have sufficient initial 
population difference s . 

During the quantum stage, the evolution is mainly due to random quantum fluctuations 
corresponding to the term 73s 2 , while coherence starts being noticeable, when the term 
72(^5 — l)w, responsible for collective effects, becomes of the same order as the quantum 
term. That is, the coherence time t co h can be estimated from the equality 

l2 (gs - l)w = 73s 2 it = t^h) , (82) 

when the quantum and collective terms coincide. This equality can hold only when gs > 1. 
Since s does not vary much during the quantum stage, the condition for the existence of 
the coherence time can be written as 

gs > 1 , (83) 
which implies that s must be positive. Equation (82) gives the coherence time 

tcoh = 7 77 ; 7 TT • (84) 

72(^0 - l)so + 73S0 + 71 (so - C) 
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In the standard situation, when 71 <C 72 and 72 ~ 73, the coherence time (84) reduces to 

1/2 

tcoh = 7 7T— • (85) 

72(330 - 1) +73 
For sufficiently strong coupling, the coherence time is 

T 2 



tcoh 



(gs > 1) • 



(86) 



To estimate the coherence time, let us take gso ~ 10. With the dephasing time T2 ~ 
10~ 13 s — 10~ 12 s, we get t co h ~ 10~ 14 s — 10~ 13 s. At the end of the quantum stage, solutions 
(79) and (80) are well approximated by the forms 



w(t coh ) ~ 2-f 3 t coh sl , s(t coh ) 



so 



(87) 



The coherence function w here is yet very small, being of order l/gs® <C 1, and the 
population difference is yet close to the initial value s - 

C. Coherent Stage 

After the coherence time, collective effects become dominant. Coherence can last during 
the time interval 

t coh <t<T 2 . (88) 
At this stage, taking into account that T 2 C^, hence, 71 <C 72, Eqs. (70) take the form 



dw . . ds 

— = -272(1 - gs)w , — = -g-f 2 w . 

These equations enjoy the exact solutions describing the superradiant pulse 

2 



w 



sech 2 

gi2/ 



t-tr 



g 



1p_ 

gi2 



tanh 



t-tr 



(89) 



(90) 



where the integration constants t p = 1/7 P and to are defined by the initial conditions (87). 
The pulse width 7 P is given by the relations 



7p = I] + ^{gi2flztcohsl , 7fl = (gS ~ 1)72 , 

which, keeping in mind that ^ztcoh «C 1, yields the pulse time 



T 2 



gso 



(gso - i) s 



1 - 



The second integration constant is the delay time 



to — tcoh + In 



7 P + 7 9 



lp-lg 



(91) 



(92) 



(93) 



corresponding to the maximum of the pulse. In view of the inequality ^tcoh *C 1, the pulse 
width can be represented as 



7 P = (gso - 1)72 + 



g 2 i2i3t coh s 2 
gso - 1 



(94) 
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where condition (83) is taken into account. Then the delay time (93) is 



to — tcoh + 




%£o ~ l) 2 

g 2 izt coh sl 



(95) 



In the case of strong coupling gs ^> 1, we have 



to — t C oh + t coh In 



2 



(96) 



irfcoh 



with r p ~ 2t coh , and the coherence time is given by Eq. (86). For 73 ~ 72 ~ 10 12 Hz — 
10 13 Hz and t coh ~ 10~ 14 s - l(r 13 s, we have ^t coh ~ 0.01 - 0.1. Then t ~ 5t coh . For 
the coupling gs 3> 1, the pulse width, as follows from Eq. (92), is r ~ T 2 /gs , which 
is inversely proportional to the dot density p, that is, it is inversely proportional to the 
number of dots taking part in the radiation process. This is a typical feature of superradiant 
emission. 

The described coherent radiation arises as a self-organized process caused by the dot 
interactions through the common radiation field. At the initial time no coherence has 
been imposed on the system, so that w(0) = 0. The radiation process is triggered by 
the transition dipolar waves, and coherence develops from the initially incoherent chaotic 
stages. This process of coherence, self-consistently arising from chaos, is the most interest- 
ing and the most difficult for description. The appearing superradiant emission is called 
pure superradiance. 

The situation is much simpler, when coherence is imposed on the system from the very 
beginning, by means of an external field, such that w(0) 7^ 0. If the system is coherent 
starting from t = 0, the incoherent stages do not exist, which implies that t co h is zero. The 
resulting coherent emission corresponds to the triggered superradiance. The superradiant 
pulse is described by the solutions of the same form (90), but with the pulse width given 
by the expression 

where w = w(0) 7^ 0. Then the pulse time is 

1 



In the case of the triggered superradiance, the pulse time depends both on the initial 
population inversion as well as on the level of the imposed coherence. 

D. Relaxation Stage 

In the time interval 




For sufficiently strong coupling, such that gs ^> 1, the pulse time becomes 



g^sl + wo 



T 2 < t < Ti , 



(97) 
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when also t 3> to, the coherent solutions (90) decay as 

(2 lp \ 2 ( 2t\ 72 ~7p , 2 7p / 2t\ 
w ~ — - exp , s ~ H exp . (98) 

\9l2j V W 912 912 V W 

Coherence dies out and the population difference relaxes to 

s ^ it » t ) , (99) 

912 

corresponding to s inverted as compared to its initial value s . For strong coupling gs ^> 1, 
when t p ~ 2t coh , t co h — l/(27 2 gso)> and 7 P ~ l29So, expression (99) equals — s , which 
implies practically complete inversion. 

E. Stationary Stage 

In the present subsection, we assume that either there is no permanent pumping, so that 
( = — 1 or that the external pumping is weak, such that \g(\ <^ 1. Then, at asymptotically 
large time 

* > , (100) 

the system tends to its stationary state. All terms of Eq. (70) play role at this stage. That 
is, the evolution is described by the equations 

dlU . . o ds , , 

= -272(1 - gs)w + 273s , — = -gi 2 w -73s - 7i(s - C) . (101) 

In order to find out the stable stationary solutions, it is necessary to resort to the Lyapunov 
stability analysis. For this purpose, we calculate the Jacobian matrix J(t) = [Jij(t)\, whose 
elements are 

d ( dw\ d ( dw\ 

Jll = d^{-dt) = 2l2{9S ~ l) ' Jl2 " Ts [it ) = 2l29W + 4735 ' 

d ( ds\ T d ( ds\ 

J21 " to, UJ = ~ 912 ' 22 " Ts UJ = " 7l ~ 72 ' 

The stationary solutions are given by the zeros of the right-hand sides of Eqs. (101). 
Evaluating the Jacobian matrix at the fixed points, we analyze the stability of the latter. 
Below, only the stable stationary solutions are presented. 

In the case when g( <^ —1, the stable stationary solutions are 

w* ~ ^ , s* ~ C f 1 T 1 ^ , (102) 

72^1 V 71I0CI/ 

which correspond to a stable node, since the eigenvalues of the Jacobian matrix, defining 
the characteristic exponents, are all negative, 

Ji - -7i - 7s , J2 ^ -272KI • 
The coherence function w is very small. 
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For weak pumping, such that \g(\ <C 1, the stable stationary solutions are 



w ~ 



7iC V 73 



7i + 73 / 7i 
7iC 



s ~ 



7i + 73 



1 - 



1 71(71 -i3)g( 

{li + Is) 2 
7i73fi'C 



(7i + 7s) 2 



(103) 



which also correspond to a stable node, as far as the characteristic exponents are 

J\ ^ -7i - 73 , ^2 - -272 • 

Because of the relations 7i 72 ~ 73, the level of coherence is very small, that is, tu* C 1. 
At the stationary stage, when there is no external pumping, coherence is practically absent, 
since it has been died out yet during the relaxation stage. The population difference at 
this stage is close to 71C/73. 

F. Pulsing Superradiance 

In the case when there is a sufficiently strong external pumping, such that g( ^> 1, the 
fixed points 

w* ~ — , s* ~ - ( 1 - , (104) 

72# g V 7i#C/ 

represent a stable focus, with the characteristic exponents 

Jl,2 - - ^ (71 + 73) ± iUeff , 

in which the effective asymptotic frequency is 

U e ff = ^2^7172 • 

The effective asymptotic frequency ou e ff defines the effective asymptotic period 

2tt 




Teff = — = 7T A /-i^. (105) 

In this regime, there occurs a series of superradiant pulses, bursting in the intervals of time 
close to the effective period (105). The total number of such pulses is of order Ti/T e //. 
In the presence of a permanent nonresonant pumping, guaranteeing the value ( ~ 1, the 
quantity r y 1 acquires the meaning of the pumping rate, which can be made comparable with 
72. Therefore, the number of pulses is of order \fg(,. Hence, there can be produced several, 
around 10, pulses. The interval between the superradiant pulses is of order T e // ~ 10~ 13 s. 

G. Numerical Illustration 



In order to illustrate the dynamics of radiation in graphical form, we calculate numer- 
ically the quantities w and s as functions of time for some parameters typical of quantum 
dots. To evaluate these parameters, we use the values from Sec. II, from where we have 
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70/72 ~ 10~ 2 — 10~ 3 and 71/72 ~ 10~ 4 — 10~ 3 . The value of the coupling parameter g, 
given in Eq. (73), is g ~ (7o/72)pA 3 . The coherence factor pA 3 ~ 10 — 10 5 . Thence, 
g ~ 10~ 2 — 10 3 . According to condition (83), well developed coherence appears when 
gs > 1. Since the initial condition s cannot be larger than 1, it should be that g > 1. 
And, as is mentioned above, for semiconductors, 73 > 72- 

Figures 1, 2, and 3 show the evolution of w — w(t) and s = s(t) as functions of 
dimensionless time t, measured in units of T2 = I/72, for several typical cases. We assume 
that at the initial time, the system is inverted, but coherence is absent and develops in a 
self-organized way. Recall that the standard mean-field, or semiclassical, approximation 
cannot describe such a self-organized appearance of coherence. Figures 1 and 2 correspond 
to the case of no external pumping, when ( = —1. The difference between these figures 
is in the value of the dynamical attenuation rate. As we see, the larger 73 decreases the 
delay time and makes the superradiant pulse strongly asymmetric. Recall that for atoms 
in free space, 73 is usually much smaller than 72, because of which superradiant pulses 
produced by atoms are more symmetric. The essential asymmetry of superradiant pulses 
is the feature typical of quantum-dot radiation. Another typical feature of the quantum- 
dot dynamics, also caused by the large rate 73, is the much faster, than for atoms in free 
space, tendency of the population difference to the stationary state. Figure 3 demonstrates 
the radiation dynamics in the case of an external pumping, when ( = 1, and there appear 
several superradiant pulses with decaying amplitude. 

10 Conclusion 

The theory of quantum-dot radiation is developed being based on microscopic equations. 
The possibility of realizing the superradiant regime is analyzed. The temporal evolution 
during all radiation stages is studied in detail. A special attention is payed to the process 
when coherence arises from an initially incoherent state. The description of this process 
is impossible by means of the standard semiclassical equations, because of which a more 
accurate method has been used in the paper, employing the stochastic mean-field approx- 
imation that has been developed earlier and applied for describing the dynamics of spin 
assemblies [5-11], Bose systems in random fields [38-40], and atomic squeezing [50]. 

It is necessary to emphasize that the radiation dynamics of quantum dots has several 
specific features distinguishing this dynamics from atomic radiation. This is connected, 
first of all, with rather different values of physical dot parameters, as compared to atomic 
parameters. Because of this, despite many analogies, the theory of dot radiation requires 
a separate investigation. The principal theoretical points that have been suggested in the 
present paper for the adequate description of dot radiation are as follows. 

(i) Because of essential current fluctuations in semiconductor, the standard semiclassical 
approximation, often used for atoms in free space, is not applicable for quantum dots. 
For the latter more elaborate techniques are required, such as the stochastic mean-field 
approximation. 

(ii) Taking into account the fluctuation of current makes the dynamic attenuation 
parameter 73 of the order or larger than 72. This is contrary to the case of atoms in 
free space, where usually 72 is the largest relaxation parameter. 

(iii) For the correct description and principal understanding of the mechanism, trig- 
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gering the beginning of the radiation process, it is important to stress the existence of 
triggering dipolar waves. 

(iv) The single-mode picture is not applicable for quantum dots. It is necessary to 
consider a bunch of transverse modes forming spatial filaments. To reduce the consideration 
to a treatable problem, it is necessary to involve some tricks, like the transverse-mode 
expansion. 

(v) The overall dynamics of dot radiation consists of several stages, which have been 
thoroughly studied and described, both analytically and numerically, for the parameters 
typical of quantum dots. 

In the dynamics of dot radiation, it is possible to distinguish the following qualitatively 
different stages. The first is the fluctuating stage lasting during the time interval < 
t < ti n t, when the radiation process is triggered by dipolar waves. At this stage, there 
is no yet sufficiently strong interaction between dots. The interaction time is of order 
t int ~ l(T 15 s - l(T 14 s. 

The second is the quantum stage in the temporal interval t int < t < t coh , when the dot 
interactions through photon exchange start playing noticeable role, but coherence has not 
yet been developed. The coherence time, required for the appearance of well developed 
coherence, is of order t co h ~ 10~ 14 s — 10 _13 s. 

Then the coherent stage comes into play in the interval t coh < t < T 2 , when the dots 
emit a coherent superradiant pulse. For the quantum dot materials, the dephasing time 
is of order T 2 ~ lCT 13 s — lCT 12 s. The maximum of the pulse occurs at the delay time 
to ~ 5t co h and the pulse duration is t p ~ 2t co f l . The pulse duration is inversely proportional 
to the dot density, that is, inversely proportional to the number of dots, involved in the 
process of radiation, which is a typical feature of superradiance. 

After the superradiant pulse is emitted, the system relaxes to an incoherent state during 
the relaxation stage in the interval T 2 < t < The population difference reverses. For the 
system of dots in a semiconducting material, the longitudinal relaxation time is T\ ~ 10~ 9 s. 
But this is not yet the final stage of evolution. 

The stationary stage is reached for t > Ti, if there is no external permanent pumping 
or the effective dot interactions are weak, so that \g(\ <C I. Then the system tends to a 
stationary incoherent state representing a stable node. 

If the system of dots is subject to a sufficiently strong external permanent pumping, 
such that \g(\ ^> 1, the regime of pulsing superradiance occurs. Then a series of about 10 
superradiant bursts can appear, flashing in the intervals of time T e ff ~ 10 _13 s. 
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Appendix 

Here, the explanation is given of the contributions coming from the dot self-action. For 
a dot located at r = 0, the self-action vector potential is 
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with the current 

J (t - = iuo [dS + (t)e- lk " r - d*S-(t)e ikor ] Q{ct - r) , (A2) 

in which S a (t) = S a (0,t). At short distance, such that k r <C 1, one has exp(ik r) ~ 
1 + ik r. Substituting the transverse delta-function (9) into the vector potential (Al), we 
keep in mind that the averaging of the dipolar tensor over spherical angles yields zero. As 
a result, the vector potential (Al) becomes 

A self (r, t) = - kl [dS + {t) + d*S-(t)] + i [dS + (t) - d*S-(t)] . (A3) 

Averaging this expression over the radial variable between the electron wavelength A e = 
27r/mc, with m being the electron mass, and the radiation wavelength A = 2n/k , and 
taking into account that A e <C A, we get the self-action potential 

Keif(t) = \ K [dS + (t) + d*S~(t)} + ^ In [dS + (t) - d*S~(t)} . (AA) 

Let us introduce the natural width 

70 = \ \A\V„ (Ab) 



fe^ln ™). (.46) 



and the Lamb shift 
2n \ uo 

Then the terms in Eqs. (14), induced by the self-action potential (A4), are 

2k S z d ■ A self = (i8 L - l2 )S- + |^ (72 + iS L )S + (A7) 

for the first of Eqs. (14) and 

k (Std + S-d*) ■ A self = 7i Q + 5 2 ) 

for the second, where 71 = 270 and 72 = 7o- In the standard situation, one has 



(AS) 



^ « 1 , — « 1 . (A9) 

The Lamd shift, without the loss of generality, can be included in the definition of the 
transition frequency ojq. 
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Figure Captions 



Fig. 1. The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T 2 ) for the attenuation parameters 
71 = 0.003,73 = 1 (measured in units of 72), for the coupling parameter g = 10, with the 
initial conditions w = 0, s = 1. 



Fig. 2. The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T 2 ) for the attenuation parameters 
71 = 0.003,73 = 10 (measured in units of 72), for the coupling parameter g = 10, with the 
initial conditions Wq = 0, Sq = 1. The larger dynamic attenuation 73 makes the pulse more 
asymmetric. 



Fig. 3. The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T2) in the case of an external 
pumping, for the parameters 71 = 10,73 = 1 (measured in units of 72), for the coupling 
parameter g = 100, with the initial conditions Wq = 0, So — 1- The coherence intensity, as 
well as population difference, exhibit five pulses with decaying amplitude. 
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Figure 1: The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T 2 ) for the attenuation parameters 
7i = 0.003,73 = 1 (measured in units of 72), for the coupling parameter g = 10, with the 
initial conditions w — 0, sq — 1. 
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Figure 2: The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T2) for the attenuation parameters 
71 = 0.003,73 = 10 (measured in units of 72), for the coupling parameter g = 10, with the 
initial conditions wq = 0, sq = 1. The larger dynamic attenuation 73 makes the pulse more 
asymmetric. 
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Figure 3: The coherence intensity w (solid line) and population difference s (dashed line) 
as functions of dimensionless time (measured in units of T 2 ) in the case of an external 
pumping, for the parameters 71 = 10, 73 = 1 (measured in units of 72), for the coupling 
parameter g = 100, with the initial conditions Wq — 0, Sq — 1. The coherence intensity, as 
well as population difference, exhibit five pulses with decaying amplitude. 
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